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Abstract. We study the effect of shear flow on the entropic Helfrich interaction in lyotropic surfactant 
smectic fluids. Arguing that flow induces an effective anisotropic surface tension in bilayers due to a 
combination of intermonolayer friction, bilayer collisions and convection, we calculate the reduction in 
fluctuations and hence the renormalised change in effective compression modulus and steady-state layer 
spacing. We demonstrate that non-permeable or slowly permeating membranes can be susceptible to a 
undulatory instability of the Helfrich-Hurault type, and speculate that such an instability could be one 
source of a transition to multilamellar vesicles. 

PACS. 64.70.Md Transitions in liquid crystals - 83.80.Qr Rheology of surfactant and micellar systems, 
associated polymers - 87.16.Dg Membranes, bilayers, and vesicles 



1 Introduction and Overview 

Lamellar systems display a wide variety of behaviours 
in the presence of shear flow. Weak flows typically align 
lamellae so that they slide over one another with layer 
normals parallel t o th e flow-gradient (Vv) direction (c 
orientation in Fig. 



1.1 



Lamellar block-copolymer melts 
undergo a series of reorientation transitions as a function 
of strain rate, frequency, and amplitude, between the c and 
a (layer normals parallel to the vorticity (uj) direction) ori- 
entations 3|. Thermotropic smectic liquid crystals 
also undergo a variety of reorientational transitions as a 
function of strain rate Lyotropic smectics, such as 
surfactant lamellar systems, can form multi-lamellar vesi- 
cles, or "onions", either above a critical strain rate (see 
for studies of SDS surfactant systems), or for very 
small (essentially zero) strain rates (see [[l0| and jyl for 
studies of AOT and ionic systems, respectively). At higher 
strain rates these onions can undergo further transitions 
to onion crystals of different s ym metries and sizes |l2|,[L3[ 
as well as to oscillatory states H] . Scattering experiments 
on some lyotropic and thermotropic smectics (main chain 
and side chain) suggest other multilamellar curved struc- 
tures, such as concentric cylinders, or "leeks" [ p^|[l^ ,[l7|. 
Onion formation and layer rearrangement have been cor- 
related with screw dislocations in DMPC/C12-E5 solutions 
]l8t . Onions have also been seen in triblock copolymer so- 
lutions n). 

The relative and absolute stability of the a and c orien- 
tations has been understood from several points of view, 
but there is no clear picture yet that indicates which mech- 
anisms operate in which systems. Cates and Milner pre- 
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Fig. 1.1. Orientations of a Lamellar Phase in Shear Flow 



dieted that shear can stabilise a layered system near the 
equilibrium sponge-lamellar phase transition pifl; s uch a 
transition has been reported experimentally ]22|,|23|. Ra- 
maswamy p4}| argued that shear flow can suppress undu- 
lations in the a orientation, leading to a collapse of the 
lamellar phase. The dependence of the critical shear rate, 
7 C on layer spacing, d and membrane viscosity rj was cal- 
culated to be 7 C ~ (kBT) 3 /r]K 2 d 3 . Such a collapse was 
seen in the flow of a thin lamellar phase of the ionic sur- 
factant AOT in brine, in a specially designed Couette cell 
assumed to shear the layers in the a orientation ]2q| . The 
dependence on layer spacing matched Ramaswamy's pre- 
diction, if the measured d-dependence of the zero shear 
viscosity (77 ~ l/<^ 2 ) was incorporated. However, viscos- 
ity measurements were performed in a standard Couette 
cell rheometer; there is some evidence that AOT orders 
in the c-orientation in such a geometry |2q| . In addition a 
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slightly different dependence of the critical shear rate on 
layer spacing was found for a lamellar phase consisting of 
an anionic surfactant, a deviation that might be explained 
by the presence of the more dominant Helfrich interaction. 
To further understand these experiments, in addition to 
disentangling the viscosity issues, it may be important to 
consider the propagation of the air/fluid interfacial surface 
tension into the bulk. 

Milner and Goulian p^jj examined the instability of 
thermotropic smectics in the c orientation, and noted that 
convection of bending fluctuations can induce compression 
and destabilises the c orientation with respect to the a ori- 
entation. Fredrickson showed that non-linearities in the 
diblock free energy can destabilise the c orientation with 
respect to the a orientation |j| . Another mechanism is vis- 
cous contrast the c phase allows the strain field to 
be concentrated in the less viscous material, which would 
be favoured at higher frequencies or strain rates. 

Although the appearance of the onion phase in shear 
flow is well-documented, there is still no reliable theoreti- 
cal framework for this phenomenon. An appealing device 
for the onion instability is an analogy of the Helfrich- 
Hurault mechanism |2^,[3^], in which an applied dilational 
strain parallel to the layer normals may be relieved by 
buckling to retain the equilibrium layer spacing. Oswald 
and Kleman [ f3l| proposed that a smectic flowing between 
misaligned plates generates defects that induce the neces- 
sary dilational strain, and Roux and co-workers used this 
argument to rationalise their results 0. The shear rate 
7 C for buckling was estimated to vary inversely with the 
square of the gap size between the plates. However, not 
onlydo experiments show that the gap size does not affect 
7 C fig] , but the predicted value for 7 C differs from the ob- 
served value by a few orders of magnitude. Wunenburger 
and co-workers Q extended this idea to consider the 
higher shear-rate transition from onions to another well- 
oriented lamellar phase. They demonstrated that, at high 
enough strain rates, an undulation could persist in the 
velocity gradient direction while a restabilisation would 
occur in the velocity direction. Reasonable comparisons 
were made between the resulting theoretical stability dia- 
grams and the experimentally observed transitions. 

In the case of thermotropic smectics, Auernhammer 
and co-workers |33) suggested that flow can induce a layer 
tilt due to brushing of molecules in adjacent layers aligned 
parallel to the normal. This induces an effective strain due 
to the layer shrinkage in the normal direction (if there 
are few defects and the necessary change in layer number 
cannot occur) that can, in principle, be large enough to 
destabilise the c orientation. Unpublished simulations [Q 
found undulations in the vorticity direction. 

In this work we consider the entropic, collision dom- 
inated membranes first studied by Helfrich |?5| , and ne- 
glect the effect of electrostatic forces. Shear flow has a 
dramatic effect on the fluctuation spectrum, and a rea- 
sonable proposition is that flow suppresses fluctuations in 
the c orientation, hence reducing intermembrane interac- 
tions and repulsion between the layers. For systems with 
many defects, this reduced repulsion would be expected 



to lead, after an appropriate time, to a reduction in mean 
layer spacin g. In deed not only was this seen by Yamamoto 
and Tanaka |36j , but also the shape of the x-ray scattering 
Bragg peak became broader and less intense. For systems 
with few defects or very slow permeation, the reduced fluc- 
tuation spectrum would lead to an effective imposed strain 
along the layer normals, which if sufficiently large could 
be relieved by buckling. This, then, could be a mechanism 
leading to onion formation, in the event that other consid- 
erations allow for onions. For example, the Gaussian cur- 
vature modulus k should not be too positive, which would 
tend to suppress structures with positive mean curvature, 
and the a orientation should also be dynamically unstable. 

Our task, then, is to calculate the suppression of layer 
undulations in the c orientation, and hence the renor- 
malised "interlayer potential" in steady state shear flow. 
In lieu of potentials, the proper starting point is the full 
dynamics of a lyotropic smectic in the presence of an 
imposed boundary stress or strain rate. Rather than tak- 
ing this approach, we assume a coarse-grained free energy, 
following Bruinsma and Rabin , as a first step towards 
understanding the effect of flow. In common with the re- 
cent work of Zilman and Granek |^H , we assume that 
bilayers oriented in the c orientation experience a flow- 
induced effective tension, although the tension that we 
propose is anisotropic and has a different sign and physical 
interpretation (Appendix ^). In Section gwe explore the 
effects of this tension on the steady-state layer spacing d 
and the compression modulus B of an equilibrium Helfrich 
stabilised lamellar phase. In Section ^ we examine these 
results in the context of an imposed shear flow. After esti- 
mating possible sources of tension we discuss the predicted 
change in layer spacing and the possibility for an instabil- 
ity, taking into account permeation and defect-creation. 
We finish with a summary in Section ^ and outline future 
work that will quantitatively address the dynamic origin 
and consequences of this tension. 

There have been several discussions of the effect of flow 
on the layer spacing in block copolymer lamellar systems, 
which focus on the change in polymer conformations due 
to the stretching imposed by flow in melts and solutions 
[ [42"|| , and layer changes have been reported in melts [ffBf . 
Our mechanism is unrelated to this one, but may be re- 
lated to recent experiments on block copolymer solutions 
that have demonstrated a change in layer spacing Jl6[ and 
onion formation as a function of flow p9(| , if those systems 
are Helfrich stabilised. 



2 The Lamellar Phase under Tension 

2.1 The Equilibrium zero-tension Lamellar Phase 

The equilibrium free energy F of a three-dimensional ly- 
otropic smectic fluctuating at fixed chemical potential pe- 
nalises gradients in the mean layer displacement u(r): 

F = i fd 3 r [K {V\uf + B{d z u) 2 ] . (2.1) 
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The elastic constants in Eq. (2.1) are the bending mod- 
ulus K and the compression modulus at fixed chemical 
potential B. For membranes that interact via collisions 
the origin of B is the steric repulsive interaction the 
collisions induce an entropic confinement pressure p, 



where a is a constant of proportionality. The calculation 
of the mean square fluctuations is straightforward: 



(* 2 (r))=X>c 



ksT 2 
4tt 3 k p ' 



(2.7) 



k B T 
L 2 p d' 



(2.2) wher e L v /a 3> 1. The constraint on the height fluctuations 
Eq. (2.6) leads to an expression for the collision length, 



where L p is the characteristic distance between collisions 
(the "patch" or "collision" length) and d is the layer spac- 
ing, from which the compression modulus may be esti- 
mated by B ~ —ddp/dd 

The patch length is the implicit lower cut-off of the 
integral in Eq. (^Tj) and is that transverse distance over 
which membranes wander before colliding, found by com- 
puting the mean fluctuations of a free membrane, governed 
by the single membrane free energy /. In the Monge gauge, 



f = \nf d 2 r(ViMr)) 2 = i^ 2 ^IVlM,(2.3) 



where h(x, y) = h(r) is local height of a given membrane, 
rather than the mean position u(r) that we introduced 
above for the coarse-grained three-dimensional free en- 
ergy. A± is the projected membrane area orthogonal to 
the mean membrane normal. L is the system size and we 
have defined the Fourier Transform h(r) — Ylg ± 
Assuming periodic boundary conditions, 



(2.4) 



where a is a molecular size and wavelengths longer than 
L p are suppressed by steric hindrance with neighbouring 
membranes. Note that 



u{x, y,z = nd) 



A(L P ) 



[h(x' - x, y' - y) - nd] 



72 ^1 



A(L P ) 
(2.5) 



ing L p to d 



The constraints of neighbouring mem- 



branes limits the mean square fluctuations to the mean 
layer spacing d 2 , 



{h 2 (r))=ad 2 



(2.6) 



L p -cd Xl r 



(2.8) 



If a = 1/3-7T 2 , as inferred by Helfrich |3f| , then c = ^At:/3. 
There is some debate about the exact value of c; for exam- 
ple, Golubovic and Lubensky |5lJ calculated c = -\/32/37r 



Combining the collision length with the pressure (Eq. 2.2 ) 
yields a compression modulus that has been calculated 
precisely as [B5|B3,M] 



B 



6S n n (k B T) 2 
n + 1 nd? 



(2.9) 



where n is the number of membranes. For n — oo Helfrich 
determined <5oo = 37r 2 /128 ~ 0.23 [p5| . Strong coupling 
perturbation calculations |32j suggest that Soc ~ 0.1 and 
for a single membrane, 5i = 7r 2 /128 ~ 0.07; both results 
are consistent with those given by Monte-Carlo simula- 
tions m . 

For flat membranes in a lamellar phase the volume 
fraction is given by <j) — t/d, where t is the layer thickness. 
For fluctuating lamellar phases of uniform thickness t, the 
mean volume fraction is 



A\ t 
d/Al 



(2.10) 



where A is total membrane area. To find the relation be- 
tween layer spacing and concentration, we can expand (for 
given mean spacing d) A ~ A±(l + ^((V ±h) 2 ) , to find 



is the average of the microscopic layer displacement over 
a pa tch area A(L p ) centred at (x,y,z = nd), so that 
Eq. (2.1) encompasses Eq. (2.3) subject to a hard wall 
constraint and K = nj d. 

We neglect the Gaussian curvature, penalised by the 
modulus S, which is valid for fixed topology. However, 
transitions to a different topology such as onions, or changes 
in the distribution of dislocations or pores |Q, may in- 
volve k. Quantities known to alter R, such as salinity, tem- 
perature, or cosurfactant, seem to influence the shear rate 
and/or the mechanism for onion formation [^fH|4^ ,^9|,^o|. 

The compression modulus may be estimated by relat- 



[1 + ±<(V^) 2 )] 

, fcsT\ / d 
Attk \ a 



(2.11a) 
(2.11b) 



Membranes dominated by electrostatic interactions are es- 
sentially flat, in which case d<$> is independent of the vol- 
ume fraction d>. Hence, deviation from this 'ideal' dilu- 



tion law according to Eq. (|2.1l[ ) suggests that the mem- 
branes are stabilised by the Helfrich interaction j54j . Use 
of Eq. (2.11), given d and t, allows k and a to be extracted 
from swelling experiments. Helfrich systems include var- 
ious dilute non-ionic and screened surfactant solutions. 
Whether or not certain dilute copolymer solutions that 
form onions under shear jl9| are stabilised by the Helfrich 
interaction is unclear. Further evidence of Helfrich stabili- 
sation is deduced from the power law behaviour exhibited 
by the x-ray structure factor |55|| , from which the bending 
modulus can be extracted. 
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2.2 Lamellar Phase subjected to a Finite Anisotropic 
Tension 

2.2.1 Equilibrium Layer Spacing 

When a Helfrich-stabilised lamellar phase is sheared in 
the c orientation, layers are convected and stretched by 
the flow, adjacent fluctuating layers collide, and flow in- 
duces the leaves of the bilayer to slide over one another. A 
rigorous study of the effect of flow on membranes should 
consider appropriate equations of motion and treat the 
dynamics explicitly. For example, the best approach, al- 
beit probably intractable, would be to coarse-grain a dy- 
namical description of coupled monolayers such as that 
of Ref. J56| to an effective two-component smectic theory, 
incorporating both flow and the normal and tangential 
collision forces. We adopt instead a 'quasi-equilibrium' ap- 
proach to gain insight into the stability of the steady states 
for membranes in flow. We propose that the resulting in- 
terplay of hydrodynamic interactions and interbilayer fric- 
tion can be described by a flow-induced effective tension 
acting parallel to the layers. 

In the presence of a tension, Seifert's (H^j self-consistent 
calculation of the potential between a wall and a mem- 
brane of a vesicle under tension might be applied to find 
the preferred separation. However, since repulsive interac- 
tions from neighbouring membranes dominate the inter- 
lamellar potential far from the unbinding transition |59| , a 
similar derivation of the layer spacing cannot be achieved. 
Thus we are led to deduce a new layer spacing without re- 
gard for the change in potential by generalising Helf rich's 
geometric argument leading to the dilution law Eq. ( 2.11 ) 
to include a lateral tension |6(J . 

To retain the asymmetry between flow and vorticity 
directions in the plane, we model the flow as an effective 
anisotropic tension applied in the flow direction. The free 
energy of a membrane subject to an anisotropic tension a 
is 



n{V\hf + a{V x hy 



(2.12) 



where a penalises short wavelength modes and the in- 
creased area due to stretching the membrane in the re- 
direction. This decreases the excess area and suppresses 
fluctuations: 




(2.13a) 



(2.13b) 



for Lp/a 3> 1. Here we assume an isotropic change in 
patch length, rather than the anisotropic renormalisation 
that would be expected; this should suffice for qualitative 
behaviour, as we show below. Assumin g th e same hard 
wall constraint as without tension, Eq. (^_^), we can cal- 
culate a renormalised mean patch length, 

rrd 2 



By the equipartition theorem, 



<(V_LMr)r> = 
k B T 



2-7TK 



ln^ 

a 



(2tt) 2 7 7r/Lp Kq\ + aql 
1- 



ln 



1 



KTT 2 



1 



1 



KTT 2 



(2.15a) 
(2.15b) 



The simultaneous condition s relat ing the patch size to the 
layer spacing via (h 2 ) (Eq. |2.13a ), and the concen tration 
to the degree of crumpling ((\7±h) 2 ) (Eq. [2.11a ) deter- 
mine the layer spacing as a function of concentration and 
tension, d((j),a), as well as the a ssociated collision length 
L p ((p, a). Hence, using Eq. (2.14) to eliminate L p ^, we find 



d = d + — — In 

4nk 



fl 



7 \ 



V 



2W1 



*(*)' 



(2.16) 



where d = d<f>/t, t = t/a, k = n/ksT and a = A^aad\jkBT 
are non-dimensionalised parameters and do is the layer 
spacing at zero tension. As expected, an increasing ten- 
sion reduces the excess area. 

For aCl we can calculate the change in the preferred 
layer spacing, 



k R T 



(i 



lT f_2_V 

i \cdo J 



47TK 4> J 



ira 



(2.17) 



Note that the linear correction is half that of a membrane 
subjected to isotropic tension. The asymptotic expression 
for high tension is 



t 



ksT c 



4ttk a V 4iraa 



(2.18) 



The relat ive reduction in la yer spacing, calculated from 
Eq. (2.16), is shown in Fig. 2.1. This applies directly to 
equilibrium systems with a mechanically- applied tension, 
such as few-layered vesicles or bilayers tethered to sur- 
faces. Tension should have the greatest effect when acting 
on the largest excess surface area, and hence we expect 
a greater decrease in layer spacing for a lower bending 
modulus and large crumpling fraction. 

As mentioned above, we have assumed that the patch 
length increases isotropically with increasing tension, and 
hence imposed a cylindrical geometry with a change in 
the lower radius cutoff tt/L p . Physically, one would ex- 
pect, at least for small tensions, that the patch size would 
actually increase anisotropically, with fluctuations in the 
transverse direction less affected by stretching. In this case 

* For an isotropic tension the swelling law is d = do + 
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Fig. 2.1. Fractional change in layer spacing d/do for the 

isotr opic (cylindrical] patch length renormalisation, from 

Eq. (2.16) and Eq. (2.11) at t = 6 as a function of induced 
tension a for small values (a) and lo g a f or all values (b). 
(c) A com parison of cylindrical (Eq. 2.16) and rectangular 
(Eqs. 2.19) patch size renormalisations, for k — 0.4, <f> — 0.1. 
Note k = K/k B T and c = 



one should impose a rectangular domain of integration, 
with a cutoff w/ L py for the q y integration, which remains 
fixed for increasing a, and a cutoff 7r/L p (0, a) < tt/L pv for 
the q x integration, which decreases for increasing a. In this 
case the simultaneous set of equations that determines the 
new layer spacing d and patch size L p is: 



ad 2 



k B T 



7r / a it I a 



dq y I dq x 



1 



47T2 J •*» J -^2 + ^2)2 + ^2 
Ti/Lpy tt/ L p (4>,cr) 
it / a tt I a 

k B T 



Wx + Ty) 



n/Lpy -k / L p ((j>,(r) 



(2.19a) 



(2.19b) 



Results of this calculation for a representative parame- 
ter set are shown in Fig. 2.1c. As expected, the inclu- 
sion of an anisotropic patch size leads to a smaller layer 
shrinkage, because the transverse fluctuations remain sub- 
stantial. Note, however, that solutions do not exist above 
for a ~ 4.0. Formally, this limit corresponds to L p — > oo 
with L y finite, and describes layers that are essentially flat 
in the tension direction while crumpled in the transverse 
direction. Microscopically, membranes are essentially in- 
compressible, so this limit is unphysical and we expect 
the layer contraction to evolve towards the isotropic patch 
calculation with increasing stress. 



3 Application to Membranes in Flow 
3.1 Physical Picture 

We consider a Hel fric h-stabilised lamellar phase with layer 
spacing do (Fig. 3.1(A)). In flow, transverse membrane 
fluctuations are suppressed, leading to fewer collisions and 
hence more "space" between layers. If there is no per- 
meation or other mechanism to change the layer number 
on the experimental time sca le, t he layer spacing remains 
fixed at its initial value (Fig. 3T(B)). At sufficiently large 



strains (Section |3.3.2 ) an instability will relieve the strain 
in favour of undulations via a non-equilibrium version of 
the Helfrich-Hurault effect. This then could be a precur- 
sor to either a stable modulated phase or an instability 
to onion formation. However, significant permeation due 
to passages, trans-membrane diffusion, or smectic defects, 
can lead to a layer spacing reduction after sufficie nt t ime 

(C). 



3.1 



for the required increase of layer number, as in Fig. 
Which effect is seen depends on the rate of application of 
flow, the permeability and the defect structure, and the 
degree of crumpling. 



(A) 
7 = 




A ± 
(B) 



actual layer spacing 
preferred layer spacing 

(C) 



7 > 




Fig. 3.1. Two different scenarios for a lamellar phase (A) sub- 
jected to flow. In (C) fluctuations are reduced and the layer 
spacing reduces due to the softened repulsive potential, while 
in (B) the equilibrium layer spacing is maintained because the 
system cannot relax (and form more layers). 



3.2 Sources of Tension 

Under the assumption that flow induces an effective ten- 
sion in a lamellar phase, we apply the results of Sec- 
This is heuristic, but we believe yields quali- 



2.2.1 



tion 

tatively useful results, as employed by Zilman and Granek 
PH . Dynamically, tension corresponds to a tangential force 
resisting changes in total membrane area that induces a 
local normal force on the membrane. Hence, we contend 
that shear flow "irons out" small scale membrane wrinkles. 
Several dynamic effects can lead to an effective tension: (i) 
shear flow distorts and stretches membrane fluctuations, 

(ii) friction between the leaves of the bilayers increases the 
drag due to intra-bilayer dissipation in undulations, and 

(iii) intermembrane collisions induce a transverse mem- 
brane force. The first two effects are non-linear and yield 
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tensions that in the simplest case scales as 7 2 , while the 
latter effect is linear in 7. All of these mechanisms are 
expected to be unimportant for flat, non-Helfrich mem- 
branes with weak equilibrium undulations. 

The importance of convection and stretching can be 
estimated by balancing the restoring force on a wrinkle 
of wavenumber q with the viscous dissipation encountered 
in the flow in a thin film to find the typical undulation 
lifetime jf(|. Fluctuations are influenced when the strain 
rate is of order the inverse lifetime of the slowest relevant 
undulations, which are those of wavelength the patch size 
L p . This leads to an estimate for the critical strain rate j c 
at which undulations are appreciably suppressed by flow 



7c 



{k B Tf' 2 
rjd 3 K 3 / 2 



(3.1) 



where 77 is the solvent viscosity. 

A more careful calculation begins with, for example, 
the Langevin equation for the affine convection of a fluc- 
tuating membrane^ in the c orientation in a shear field 
v = 7ZX, 



The noise £(r, t) may or may not be related to r(r — r') 
through the fluctuation-dissipation theor em, depending 
on the strength of the applied shear. Eq. ( |3.2| ) is a varia- 
tion of an anisotropic Burger's equation with a convective 
non-linearity Q. Upon coarse-graining Eq. ( [3.2] ) up to 
the collision length L p , the convective non-linearity gen- 
erates a dynamic response analogous to an effective scale- 
dependent tension, 

ki 

- [r q±Kq i + r x (L p , »)rq 2 x ] V, (3-4) 

where r x (L p , /i) depends on the particular relaxation mech- 
anism. Thus, 



r x (L p ,n) 



(3.5) 



For example, for isolated membranes we find a "tension" 
of order 



L 5 p v 2 i 2 q±- 



(3.6) 



h(v,t) 



-K J dVr(r-r') V 4 /i(r',t) +£(r,i). 



(3.2) 



The kinetic coefficient r (r — r') depends on the details 
of the fluid-membrane coupling, with the form (in Fourier 
space) r q± ~ j]~ 1 l l/+1 q±. Here we give three examples of 
relaxation mechanisms that alter the exponent v and the 
associated length scale I: 




permeable 
isolated 

confined fluid/squeezing. 



(3.3) 



The permeation length scale £ depends on the size and 
the density of the pores and the membrane thicknes s]* * *| . 
For wavevectors q±C ^ 1 ^ q±L p the membrane may 
be considered impermeable. In the confined fluid regime 
(qj_d 1), squeezing of solvent within the confines of 
the surrounding membranes leads to v = 2 plj ]. The case 
of v = — 1 describes the hydrodynamic interaction of an 
isolated membrane in solvent {q±d <C 1) |32|, with a length 
scale fixed by the wavevector. For permeable membranes, 
the characteristic length ( < d, with relevant wavevectors 
q±e <C 1 <C q±(, where e is some small length scale at 
which this dynamical description breaks down. 

** In the dynamical description of the long wavelength vari- 
able u(r) the convective term is yzdxii, appropriate for wave- 
lengths larger than the layer spacing and for fluctuation am- 
plitudes u smaller than the layer spacing ]27|,^o|. 
* * * A simple model of flow through circular pores of width w 
leads to £ ~ tip 2 /w 4 , where ip is the mean pore separation 
within the membranes. 



This calculation, along with results for the other relax- 
ation mechanisms will be presented elsewhere |65f . A sim- 
ilar quadratic scaling was also estimated by Zilman and 
Granek based on energetic arguments. It is impor- 
tant to recognize that, although the "tension" above ap- 
plies, strictly, only to wavelengths of order the collision 
length, it is generated at all wavelengths larger than the 
smallest cutoff and grows during the coarse-graining pro- 
cedure. We have chosen to model this tension, which is 
wavenumber-dependent, as an average value that applies 
for all wavenumbers. This certainly changes any quanti- 
tative predictions, but does not influence the qualitative 
aspects of our results. 

Compariso n of the bending and tension relaxation mech- 
anisms of Eq. (3.4) leads to an estimate similar to Eq. (3.1) 
for the strain rate 7 C at which fluctuations are significantly 
suppressed 



7c 



k B T 
ijd 3 



(3.7) 



If we consider a more detailed picture of the bilayers 
there are at least two more conceivable sources of tension. 
The two leaves of the bilayer exert transverse intermono- 
layer forces as they slide over each other. For flat mem- 
branes this simply renormalises the solvent viscosity; for 
highly crumpled membranes the viscosity is significantly 
enhanced, 



Veff -V 



1 



^((V^) 2 ) 
drj 



(3.8) 



where r\h is the bulk viscosity of the bilayer material (typ- 
ically alkyl tails). In this highly crumpled case there may 
also be an effective tension, or restoring force against the 
dissipation induced by bilayers with enhanced crumpling. 
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A third contribution comes from the tangential force 
incurred in interbilayer collisions. Because the bilayers re- 
tain their integrity they are moving at the velocity of the 
mean bilayer position, rather than at the velocity that 
would be determined by affine flow. The enhanced local 
shear rate near the collision event, over a length of or- 
der a bilayer thickness rather than the interlayer spacing, 
leads to an excess tangential force per unit length due to 
collisions, or equivalcntly a tension, 



(d — t) . A con 
(Jb-b ~ V — 7 — 1—? — • 



(3.9) 



Here, A con is the collision contact area, expected to be 
quite small. 



3.3 Scenarios in Flow 



3.3.1 Layer Spacing Change 



As mentioned in Section 2.2.1, an effective tension is ex- 



pected to change the mean layer spacing in smectics with 
sufficient permeability. Specifically, either dislocation loops 
must be easily generated or the lipid permeability must 
be high enough for new layers to form before any incu- 
bation time for an instability or other transformation, 
such as to onions. Note that simple solvent permeability is 
not enough. Hence, we expect that highly defective mem- 
branes are most likely to exhibit a change in layer spacing 
under flow. An example of such a system is C12-E5, which 
is known to be highly permeable in regions of the phase 
diagram |47j] and dominated by the Helfrich interaction. 
In this case o ne ex pects a change in layer spacing gov- 
erned by Eq. (2.17), with the tension a given by a com- 
bination of collisions (Eq. 3J3) and convective effects (Eq. 
3.6), schematically of the form 



a = 777 [A + Bj^ 1 ] 



(3.10) 



where A and B are poorly known. 

Yamamoto and Tanaka found a flow-induced change 
of layer spacing in a C12-E5 lamellar system only a few de- 
grees away from the equilibrium lamellar-to-sponge tran- 
sition. This does not, of course, necessarily validate our 
mechanism, since any mechanism that predicts the gener- 
ation of new layers will be more susceptible in the presence 
of defects. In the vicinity of this temperature shear thin- 
ning behaviour was also observed, indicating that flow- 
induced change of structure was indeed occurring. They 
did not, unfortunately, report the quantitative dependence 
of layer spacing on strain rate, so it is difficult to make a 
precise comparison. Shear thinning has been seen in other 
layered systems (for example p2f ) but as far as we know 
has not been correlated with a change of layer spacing. 

The same C12-E5 system has also been observed to form 
onions under shear in other regions of the phase diagram 
The dependence of structure formation on concen- 
tration, strain rate, shear history, and time would help 
to discern whether or not the mechanism we propose is 



relevant for this system. Since flow modifies the fluctua- 
tion spectrum and the steric repulsion, the shape of the 
static structure factor should also change, which possibly 
explains the rounding of the Bragg peak observed by Ya- 
mamoto and Tanaka p6| . 

Note that, although defect generation is a necessary 
condition for changing the layer number, the presence of 
too many defects should invalidate layer suppression ef- 
fects. For example, significant porosity at a scale less than 
the patch length implies a lack of integrity of the mem- 
brane. In this case, although on symmetry grounds the 
long wavelength theory would be the standard two fluid 
smectic hydrodynamics, the approximation of the micro- 
scopic membrane theory by a simple Helfrich Hamiltonian 
would break down. For example, we would not expect the 
theory to apply to the defect-ridden SDS/decanol/water 
lamellar system achieved by low concentrations of cosur- 
factant, that was observed to undergo transitions between 
a and c orientations f67|| . In addition there have been sug- 
gestions that in some lyotropic lamellar systems the defect 
density of highly porous systems (before any discontinuous 
transition) can change continuously with the flow Jl7[ . 

3.3.2 Instability 

The alternative scenario is that in the absence of perme- 
ation, or if permeation is slow enough, the layer spacing 
cannot change from do to d. Instead however, flow would 
induce an effective dilational strain 7, given by 



7 



do — d 



(3.11) 



At sufficiently high flow, and hence strain, the system 
would become unstable to an analogue of the familiar 
Helfrich-Hurault effect, in which a smectic- A liquid crys- 
tal exchanges dilational energy for buckling energy for a 
large enough dilational strain. In that case, minimisation 
of the smectic elastic free energy determines the critical 
strain 7* at the onset of the instability and the undulation 
wavevector q* ± @|o). 

The analogous effect for the sheared smectic is that 
fluctuations are suppressed by the flow, leading to a smaller 
"preferred" layer spacing; in the absence of layer-creation 
the effect is an induced strain, without actually applying a 
dilational pressure. We can gain a heuristic understanding 
of the effect by performing an equilibrium calculation for 
the critical strain, and compare it to the induced strain. 
The free energy of the smectic- A (Eq. 2.1) is adapted |30] 
to include the non-linearity required by rotational invari- 
ance^, as well as a long wavelength symmetry breaking 
anisotropic tensi on arising from coarse graining the free 
energy (Eq. p. 12 ) up to the patch length, L p : 



F 



d A r 



B 



d z u - i(Viu) 



+K{V\uy + -(V x u) 



(3.12) 



The correct form of Eq. (3.12) involves |(Vtt) 2 rather than 



J_U) 



This shifts 7* and q* by higher orders of A/L. 
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In contrast to a positive tension acting on the ripples, im- 
posing the constraints of fixed total and projected mem- 
brane areas, as in Ref. 0], leads to a negative tension (or 
lateral compression) that causes the membrane to buckle 
under sufficient shear (see Appendix A). 

To calculate the stability, we consider a small pertur- 
bation Su around a strain induced displacement 7, 



u = 72 + Su sin q z z cos q x x cos q y y. 



(3.13) 



Upon linearising in Su, we find an undulatory instability 
with an undulation wavevector given by 



0. 



and a critical strain given by 



2Att 



7 



_ do _ 1 _ 

~ d* L 



(3.14) 



(3.15) 



where A = y K/B is the penetration length and L is the 
sample thickness. If we ignore any renormalisation of B 
due to the change in undulation spectrum, then the cri tical 
tension a c can be found, once 7* is known, from Fig. |2.l| . 
However, the suppression of fluctuations and the reduc- 
tion in collisions diminishes B (see Appendix [^), which 
increases the critical strain and hence also increases the 
critical tension. An applied tension renormalised B ac- 
cording to (Appendix M) 



B = 



9tt 2 (k B T) 2 k B Ta 
"64 Kd 3 
1 /5tt 

1671- 2 It 



7T 



Annd 
5 9 \ d 
~ ~ 16 



(3.16) 



* 2 + 



K 



I f the strain is in the linear regime, as described by 
Eq. ( 2.17 ), assuming no permeation and a <C do <C L the 
effect on the reduction of B is small and an instabil ity is 
indu c ed fo r lar ge eno ugh tension given by (from Eqs. 2.17 , 
3~Tl"l, Ol and |3T6|) 



a > a c 



128tt 2 k 2 <; 
tk B TL 



1 



tk B T 
Ann(j)do 



(3.17) 



21 J, t ~ 3nm, L ~ 1mm and 
Ak B T then a c is of the order 



If = 0.2, k B T = 4-10 
letting k range from 0.2 

10~ 8 — 10~ 5 Jm -2 . In scaled parameters, this corresponds 
to a ~ 10~ 3 — 10 _1 , and a very small critical strain isnec- 
essary, corresponding to the linear regime of Fig 



2.1. If 



<t ~ r/jdo ( owin g to the difficulty of estimating the prefac- 
tors in Eq. 3.1C ) and r\ ~ 1 mPa-s then j c ~ 10 3 — 10 5 s -1 
for same range of parameters. To achieve a value compa- 
rable with experiments, Zilman and Granek rectified 
a similar discrepancy for the critical strain rate by replac- 
ing the solvent viscosity by the measured viscosity of the 
lamellar phase so that j c ~ 1 s -1 for low k. On the other 
hand, reasonable experimental shear rates are obtained for 
small concentrations (<j> ~ 0.01). As in previous works jHJ 
[l6|,|d, 32 1 , we have taken L to be the system size. Another 



possibility is that L corresponds to the lamellar grain size, 
whose scale is set by the defect density. In this case the 
reduction in B and hence the change in critical strain and 
the strain rate will be significant. 

The observation that the layer spacing does not change 
prior to onion formation does not seem to have been con- 
firmed in experiment. However, the fact that most of the 
lamellar phases that undergo the onion transition seem 
to exhibit Newtonian or only weakly shear thinning flow 
(which is probably due to improvement in orientation of 
the layers or a reduction in defect density) indicates that 
flow is not bringing about such a structural change. 

As might be expected, the membrane is most suscep- 
tible to undulations transverse to the applied flow direc- 
tion. Hence, this instability is first to a stripe-like undula- 
tion, with stripes parallel to the flow direction, or equiv- 
alently with the wavevector in the vorticity direction. On 
the other hand, an isotropic tension that penalises Fourier 
modes in both directions favours square lattice buckling in 
preference to stripes. Although, as for an envisaged square 
lattice buckling, the mechanism from the stripe state to 
further transitions such as onions is unclear. Our result is 
consistent with recent work (unpublished) by Tanaka et 
al. that shows such an undulation in a C12-E5 surfactant 
lamellar phase |66|. Similar behaviour has recently been 
demonstrated with simulations of a thermotropic lamellar 
phase p4j . The latter does not correspond to a Helfrich- 
stabilised lamellar system, while the experiments probably 
do. Recent experiments on C10-E5 by Zipfel et al. revealed 
a kinetic intermediate between the c-oriented lamellar and 
onion phases, compatible with cylindrical multilamellar 
"leeks" " 



4.s 



Such leeks are compatible with an initial un- 
dulation in the layers, with wavevector parallel to the vor- 
ticity direction, that subsequently breaks the layers and 
stabilises in the cylindrical symmetry. We also note that 
the transient transitions to leeks and then onions occurred 
at roughly the same strain (or order a few thousand strain 
units), for different strain rates. 

We emphasise that we have ignored Gaussian curva- 
ture, which is likely to influence transitions involving layer 
topology changes under shear. In equilibrium, reduction 
of salinity (l0[^l or temperature ||| or cosurfactant [[70] 
are all known to decrease the Gaussian curvature modu- 
lus R, thus favouring phases consisting of spherical over 
bicontinuous structures. A C\q E3 /water solution devel- 
oped an onion transition upon reduction of the temper- 
ature at an imposed shear rate which is consistent 
with the known decrease of R with decreasing tempera- 
ture in this system. Similarly, a low salt ionic surfactant 
AOT lamellar system showed no critical strain rate [ |To[ 
fl8t , but rather onion formation at a characteristic critical 
applied strain with a time that increased dramatically for 
higher salt concentration and thus larger R. In this par- 
ticular case, onion formation was independent of strain 
rate, suggesting a mechanism different from that we have 
proposed. Finally, at high enough concentrations of co- 
surfactant, the SDS/decanol/water system seemed to be 
prevented from forming onions f50[ , even though the sys- 
tem was less permeable. Here we suggest the increase in 
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Gaussian curvature may also be correlated with this ob- 
servation. 



4 Summary 

We have studied the suppression of the undulations in a 
Hclfrich-stabilised lyotropic lamellar phase in shear flow, 
by crudely modelling the flow as an effective anisotropic 
tension. This decreases the intermembrane compression 
modulus B due to a reduction in fluctuations, and cor- 
respondingly changes the structure factor. There are two 
general consequences, depending on the permeability of 
the lamellar phase (under shear). 

1. If permeation or defects allow the generation of new 
layers, such a system would eventually attain a new 
layer spacing, consistent with the flow-induced reduc- 
tion of collisions. 

2. If the system cannot change the number of layers, or 
the process is very slow, then either: 

(a) The system can maintain the original layer spacing 
for low strain rates, with a concomitant effective 
strain due to the deviation from the preferred sta- 
ble steady state. 

(b) For higher strain rates, an instability to undula- 
tions along the vorticity direction can relieve the 
effective induced strain. A linear analysis does not 
allow us to determine whether or not such a state 
would result in a stable undulatory phase (as re- 
ported in simulations 34 and experiments |6(|), 
become unstable to ripping layers to produce ei- 
ther cylinders with axes parallel to the velocity di- 
rection, as seen in a number of systems, or produce 
onions. In such a situation the layer spacing prior 
to forming the new state (striped undulation, cylin- 
ders, or onions) would not be expected to show a 
change in layer spacing. 

Although the reduction of B has a very small effect 
on the critical strain (or equivalently, the critical applied 
tension or strain rate) at which an undulatory instabil- 
ity may occur, the reduction in repulsion allows any at- 
tractive potential to become more significant. For highly 
swollen lamellar phases close to an equilibrium unbinding 
transition |39) there is the possibility of a flow-induced un- 
binding, which would probably be manifested in solvent 
expulsion or macroscopic phase separation. 

In addition, we have discussed possible sources for the 
tension, including bilayer collisions, intra-bilayer dissipa- 
tion and convection. The latter mechanism leads to a scale- 
dependent "tension" upon coarse-graining from the mi- 
croscopic layer position h(x,y) to the mean smectic layer 
displacement u(x, y, z). 

One implication of our result is that, in principle, lamel- 
lae and onions in flow are generally expected to have dif- 
ferent layer spacings and hence different concentrations. 
Onions, particularly the large ones that appear at low 
shear rates and are the likely candidates for coexistence 
with lamellae [|l6| , are very close in free energy and hence 



layer spacing to the equilibrium lamellar phase, while the 
lamellae subject to shear flow could be subject to strong 
modification of their fluctuation spectrum. This would im- 
ply slow kinetics for the eventual transformation to onions, 
due to concentration changes, and non-trivial shapes for 
the measured flow curves, or 'plateaus', in regions of macro- 
scopic coexistence J72|. 

Our work has focused on an instability to onion forma- 
tion. We have not determined whether or not onions are 
metastable above a given strain rate, in which case the in- 
stability would be the analogue of an equilibrium spinodal. 
Indeed, it would be very interesting to study experimen- 
tally the history- and time-dependence upon cycling the 
strain rate or shear stress (both increasing and decreas- 
ing) through the lamellar-onion transition. Experiments 
on SDS/decanol/dodecane/ water have indicated regions 
of concentration, or layer spacing and degree of crumpling, 
in which the transition is apparently continuous (low con- 
centration or larger equilibrium layer spacing) and first 
order (high concentration and smaller layer spacing) HQ. 
In addition, the history dependence may alter the eventual 
point of instability. For example, if the strain rate is in- 
creased slowly compared to the time scale associated with 
permeation, a reduction in layer spacing may be observed 
which would yield a higher critical strain rate. Further, if 
a system displays a shear induced reduction in layer spac- 
ing, onions might be formed if a strain rate larger than 
the critical strain rate is applied faster than permeation 
effects can occur. 

We have assumed that flow acts as an effective ten- 
sion. This is obviously quite crude, and forthcoming work 
will study the dynamics of individual membranes in flow. 
For membranes in the a-orientation, as studied by Ra- 
maswamy to describe a layer collapse transition p^,p5[, 
the convective term is linear, while in the c orientation. 



the convective term is non-linear, as in Eq. (3.2). On sym- 
metry grounds, upon coarse-graining the fluctuations to 
a scale of order the patch size L p between collisions, the 
convective non-linearity generates a tension-li ke r estoring 
term in the long-wavelength dynamics, Eq. (3.4), which 
depends on the single membrane relaxation mechanism 
and is also, strictly, dependent on length scale and wavenum- 
ber in a very different form than a conventional "tension" . 

We thank H Tanaka, T Kato, J Penfold, D Roux, AJ Bray, 
TCB McLeish, D Bonn, J Leng, and K Kremer for helpful 
discussions. 



A Relation to Previous Work 

In a related work, Zilman and Granek (ZG) stud- 
ied the effect of shear flow on the c-orientation of lamel- 
lar surfactant systems, to address the lamellar-to-onion 
transition. Our treatment differs from theirs in several re- 
spects. ZG advocate the physical picture that, for cylin- 
drical Couette flow, layers retain their integrity so that 
the projected area is fixed by the experimental geome- 
try. Shear flow removes small scale fluctuations and hence 
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would stretch the projected area of a free membrane. If 
this projected area is constrained, then the membrane can 
only undergo macroscopic buckling to redistribute undu- 
lations from small scale crumpling to larger scale undula- 
tions. Hence ZG model the coarse-grained membrane as 
experiencing a negative buckling tension. Our approach is 
quite different, and consists in examining the effect of a 
positive tension-like quantity on the underlying crumpling 
spectrum and hence the steady state layer spacing. 

The ZG picture is certainly relevant if for example, 
most layers in a cylindrical Couette geometry have cylin- 
drical topology with edges only at the cylinder ends. On 
the other hand, if most layers in the system do not encircle 
the centre Couette cylinder but end in edges, then they 
need not maintain a given projected area, but can adjust 
by moving defects. We feel that this relaxed constraint 
applies more generally than the fixed projected area con- 
straint, but we acknowledge that line dislocation densities 
are notoriously difficult to determine experimentally. 

Before outlining other differences we review the differ- 
ent layer variables for a lyotropic smectic: 



h(x,y) 
u(x,y,z) 



microscopic layer position 
(a < dx, dy < L) 

smectic layer displacement 
(L p < dx, dy, dz < L), 



(A.la) 
(A.lb) 



where L is the system size, and dx, dy, and dz refer to dif- 
ferences of the independent variables x,y,z. The broken 
symmetry variable u(x,y,z) describes the average layer 
displacement. In the undcformed state a layer is taken 
to be "flat", u(x,y,z) = 0. In principle, u{x,y,z) can be 
obtained by coarse-graining the highly crumpled micro- 
scopic variable h(x,y) up to a length scale roughly of or- 
der the collision length L p ; hence u(x, y, z) is defined only 
on length scales larger than this cutoff. Conversely, h(x, y) 
is defined down to a microscopic scale a, of order a sur- 
factant head diameter. In performing the coarse-graining 
information about smaller length scales is retained, and 
resides in the other hydrodynamic variable, the concentra- 
tion, which essentially measures the degree of crumpling at 
length scales smaller than the collision length L p . Hence 
the free energy for lyotropic smectics is most naturally 
defined in terms of average layer displacement and con- 
centration changes, although other combinations of these 
basic degrees of freedom are, of course, possible. 

ZG draw an important distinction between different 
"projected areas" . ZG define the geometric area A geom as 
the conventional projected area A±_ imposed by sample ge- 
ometry. For parallel plate rheometer A± is the plate area, 
while for a Couette rheometer of cylinder height h, A± 
is the cylindrical area 2irRh at a given radius R. ZG dis- 
tinguish this area from a "physical" projected area A p h ys , 
which can be interpreted as the area of the coarse-grained 
smectic variable u; i.e. (in the Monge gauge) 



iphys 



1 + \ (V±u) 



(A.2) 



This corresponds to the "constraint" introduced in Eq. (ZG- 
19), which is in fact not a constraint but a definition, as 
follows by examining their Figure 3 in conjunction with 
their discussion introducing physical projected area (A a 
in their notation, in Figure 3 of their work). 

By modifying the derivation of the Helfrich interaction 
ZG derive an expression for the macroscopic free energy 
Fzg[ u ,A] (Eq. ZG-14a) of a lyotropic smectic as a func- 
tion of the mean layer displacement variable u{x, y, z) and 
the local change in physical projected area, 



A 



SA, 



phys 



4(0) 
phys 



(A.3) 



However, A cannot be an independent local hydrodynamic 
variable, like concentration (or volume fraction) in the 
standard two-fluid smectic desc ription, because it is rigor- 
ously defined through Eq. (A.2) as a non-local function of 
u(r). Hence, A does not contain, a priori, any information 
about small scale crumpling or concentration. 

Although ZG claim that A may be eliminated in favour 
of the local volume fraction </>, and then "integrated out" 
to recover the free energy as a function of u(r) at con- 
stant chemical potential, for which layer compressions are 
penalised by the usual modulus B (ZG:Appendix B), this 
is not cor rect. The local volume fraction is defined by 
(Eqs. [Hoi and ZG-B.l) 



tA 



(A.4) 



ZG use this relation, at fixed membrane area A rather than 
fixed projected area A± (Eq. ZG-B2), to relate changes in 
layer spacing, concentration, and A±: 



Sd 
~d 



SA ± 
' A, ' 



(A.5) 



They assert that A± may be tak en to be the physical pro- 
jected area A p h ys and apply Eq. ( |A.3| ) to relate 6<j>, Sd and 
A, but this is in consistent with the geometric definition 
above (Eqs. A.2 and ZG-19). They have already defined 
Sd in terms of u, using the geometric relation (ZG-13) 



Sd du 
d dz 



(A.6) 



corr ect t o ord er 8 2 for small rotations 9. Combining Eqs. (A.2 
A.5 and A.6) gives 



du 
dz 



(A.7) 



Thus a small rigid body rotation^ about the y axis, 
9 2 leads to the unphysical result 



u >c 
lefi 

2 C 



— XV + 2% 

Similarly, in the absence of the buckling tension a 
(which they argue vanishes for zero shear) the final free 
energy ZG use (Eq. ZG-21a), obtained from Eq. (Z G-14a ) 
by implementing the definition of physical area (Eq. A.2 ), 



In general, u = — x sin 9 + z(l — cos ( 
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is not invariant under small uniform rotations for a = 0. 
ZG also consider a scenario in which the membrane area 
can vary at fixed layer number. In this case they minimise 
Fzg[u, A] over A and use the resulting free energy to de- 
termine stability; however, this i s inc onsistent with the 
definition of A, according to Eq. (A. 2). 
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A final probl em concerns the non-linear analysis. The correct form of the rotational invariant term in the free 
energy (Eq. 3.12 and ZG-14a) involves ^(Vti) 2 rather than ^(Vj_m) 2 . The effect of this approximation on the linear 
stability is negligible; the critical strain and wavevector shift by order X/L where A is the penetration length. However, 
a substantial correction is expected for a non-linear analysis of the buckled state, as in ZG-Section 4.2. 



B Calculation of Tension-Renormalised Compression Modulus B(a) 



To compute B we follow the mean field variational approach of Rabin and Bruinsma J4(J, who calculated the free 
energy difference between a confined system and a corresponding set of free layers, and add an anisotropic tension. 
A related calculation has been performed by Lubensky et al. J73]. First, we replace the free energy of interacting 
membranes by an equivalent single membrane, 



fs(v3» + {a + ^){v x hy + rh 2 



(B.l) 



where r and \i are variational parameters that maintain the equilibrium layer spacing, while a penalises changes in 
excess area due to undulations in the ir-direction. Assuming constant membrane thickness and area per molecule, \i 
is an effective chemical potential. For convenience, we consider a reservoir that admits material along the x direction 
and also penalises excess area. The last term is a harmonic potential that mimics membrane interactions, and thus 
determines the mean layer spacing. 

The equilibrium distribution of membrane height fluctuations is given by 



P±(h) 



exp 



■i0£ q >9x+(/*M<g+r)|/ l c 



N 



where (3 = 1/fcsT and 



(B.2) 



(B.3) 



By the equipartition theorem the mean layer fluctuations are 



(h 2 (r)) 



d 2 q± 



An 2 J K q A ± +{n+a)q2+r' 



(B.4) 



The relation of (h 2 ) to layer spacing (Eq. 2.6) determines r as a function of /i, while jj, is finally determined by 
minimising the mean membrane free energy, given by (using Eq. B.2) 



F= (f)+k B T(lnp ± ) 
1 



E 



(nqi + aq 2 x )(\h (ix \ 2 } + k B T (in 



exp - §/? [nqj+jLi+p) \h q 
N 



(B.5a) 
(B.5b) 



The free energy difference between free and confined membranes is 

-W x + n , \Kqi+(fi+a) ql+r 



AF 



Kql+in+a) ql+r 



In 



nq\ + aq% 



(B.6) 



Expressing F to second order in (jj, + a) 



r(n + a) = r + r^/i + a) + -r^'( M + a) 2 + 



we find, using Eqs. (2.6, B.4), 
1 

r n = 



(3 2 K(8ad 2 Y 



r' = - 



i 



8iraf3 2 Kd 2 



r" — 
1 o — 



(3 2 K 



3 27 
16 ~ 16tt 2 



(B.7) 



(B. 
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For a — 0, Tq is the same while i~"o = — l/4 ira0 2 KcP (correcting a sign error in p0[ ) and _T"o 
Substituting Eqs. (B/7) and ( |B.8| ) into Eq. ( |B.6| ) gives, after some calculation, 



1 



1 



A± 128 nf3 2 ad 2 8\16 



1/27 l\arf : 



IQirfin 



In — -4 In 2- -it 2 ad 2 (fj, + a) 
Kl 2 



The compression modulus at constant chemical potential is given by 

MS)©- 

where Z\_F is the minimum free energy with respect to the chemical potential, yielding 



B 



9tt 2 (k B T) 2 k B To 1 
"64 nd? Annd ~ 16tt 2 



57T 

T 



16/ k C 



(B.9) 
(B.10) 
(B.H) 



As expected, Seifert's calculation of the potential between a membrane and vesicle in the presence of a small tension 
reduces to a similar result 1581 



B = 



9tt 2 (k B T) 2 irk B Ta 



64 nd 3 



32nd 



(B.12) 
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